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Construction of a Solar Power Satellite (SPS) would require the injection of 
large quantities of propellant to transport material from Low Garth Orbit (LEO) to 
the construction site at Geostationary Earth Orbit (GEO). This injection, in the 
form of ~10 32 , 2 KeV argon ions (and associated electrons) per SPS, is comparable to 
the content of the plasmasphere (~10 3 * ions). In addition to the~mass_deposited, 
this represents a considerable injection of energy. 

The injection is examined in terms of a simple model for the expansion of the 
beam plasma. General features of the subsequent magnetospheric convection of the 
argon are also examined. 


INTRODUCTION 


In recent years a large scale energy system, the Satellite Power System (SPS), 
has received considerable attention from the scientific and technical community. 

The basic concept for SPS is as follows: Large (10 km * 5 km) platforms would be 

constructed in geostationary earth orbit (GEO) to collect solar energy. This energy 
would be converted into microwaves and beamed down to Earth, received by a recti- 
fying antenna and fed into the power grid. 

In 1978, Rockwell International did a system definition study (ref. 1) in which 
a 5 GW (at Earth interface) reference system was developed. The transportation 
component would mandate the construction of several reusable heavy lift launch 
vehicles ( HLLV ) to haul material into low Earth orbit (LEO). From LEO the cargo 
would be shuttled to the construction site GEO in a fleet of electric orbit transfer 
vehicles (EOTV). 

The EOTVs would be solar powered and propelled by argon ion thrusters. Ion 
thrusters have some advantages over chemical rockets. They can deliver a sustained, 
steady thrust. Also, the ion thruster propellant velocity is much greater than for 
chemical thrusters, therefore much less mass need be injected to move an equal 
amount of cargo from LEO to GEO with ion thrusters. 


In spite of the great efficiency of ion propulsion, due to the great mass that 
must be transported, enormous quantities of energetic argon and the associated 
(thermal) electrons would be injected into the environment. Such a large scale 
Injection of plasma into the magnetosphere is likey to have a global impact on 
Earth’s magnetospheric morphology and dynamics. In addition, due to the highly ani- 
sotropic velocity distribution of the argon, this represents a considerable injec- 
tion ot free energy. Numerous processes should transfer a large portion of the 
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Injection energy to the magnetospherlc system. 

ion thruster technology Is still developing and so the parameters for the 
thruster, in question are uncertain. The Rockwell reference system described an EOTV 
propelled by an ion thruster with a grid potential limit bf 2 keV, to avoid arcing 
to the background plasma. It Would operate with an ion beam current of 1904 amps, 

have a- radius a 38 cm, and develop 69.7 Nt of thrust. Eighty such thrusters would 

equip each EOTV, in four groups of 20, with 16 active and 4 spares. 

The thruster produces two distinct plasmas: the beam plasma and the thermal 

plasma produced by charge exchange between the beam plasma and escaping un-ionlzed 
argon.- The number of charge-exchange ions produced per second is given by Kaufmann 
(ref. 2) to be: 
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where = charge exchange cross Section = 2 x 10“ 19 m 2 
Jjj = beam current ~ 1904 A 
V 0 = (SKT/nm) 1 / 2 
R = beam radius = . 38 m 
n = fraction of propellant ionized 

Carruth and Brady (ref. 3) state that in experiments with a 900-series, Hughes 
mercury ion thruster approximately 90% of the propellant is ionized. The remaining 
10% escapes through the optics in the form of neutral mercury. Therefore n is 
assumed to be 0.9 and KT = 10 eV (ref. 4), which gives N = 6.07 x 10 20 Ar + s _1 . 

This represents about 5.6% of the beam current. 


and Katz (ref. 5), and Carruth and Brady (ref. 3) report that laboratory 
tests Show the charge— exchange plasma near the thruster moves radially outward from 
the thruster beam. This thermal plasma will be injected into space with essentially 
the EOTV's orbital velocity. As in the barium release experiments (ref. 6), the 
plasma is expected to expand until B ~ 1, at which point the expansion perpendicular 
to B is stopped by the field. 

This thermal argon plasma, apart from the beam plasma, would be in itself a 
considerable addition to the thermal heavy ion population, especially in the plasma- 
sphere. The remaining unionized argon would be subject to charge exchange and pho^- 
toionizatlon, the latter of which has an e— folding production rate given by Siscoe 
and Mukherjee (ref. 7) to be 4.5 * 10~ 4 s” 1 . This allows many of the fast (charge- 
exchange) neutrals to escape, while trapping the thermal neutrals in, oi near, the 
plasmasphere. 


TABLE I. - SPS PARAMETERS 


SPS mass ~ 5 * 10 7 Kg 
EOTV — - LEO departure ~ 6.7 * 10^ Kg 
Cargo ~ 5 x lo & Kg 

Propellant 5.5 x io 5 Kg 

LEO ♦ GEO + LEO trip time = 130 days 


The quantities of mass Involved In the reference system are given In table l. 

Tn hulld a 5 GW station of a mass * 5 x 10? Kg, one needs * 10 EOTV flights which 
would Inject 5 * 10 6 Kg Of 2 keV Ar + and ~5 x 10 5 Kg of thermal- Ar (along with t e 
associated electrons) Into the magnetosphere. ^|t two staUons are built- 

ner vear this gives an average Injection rate of 5.3 x ip 40 Ar s ♦ mis is com 
Sriu to Ehe polar wind Infection rate of -3 . 10“ s' and equal to the plasna 
sheet loss rate (ref. 8>- The average rate of energy Injection (In the for. °{ 

2 keV Ar + ) is ~10 9 watts, while an average substorra deposites 10 10 watts into 

the ionosphere (ref. 9). This energy will, however, he distributed over a seller 
area and so power densities could be similar to auroral power densities. 

Needless to say, the scope of the questions invoLved in such an injection is 
extensive. S This ££« will concern Itself mainly with two topics: the Injection of 

the energetic Ions (beam plasma dynamics) and the subsequent convection of the beam 
ions in the magnetosphere. 


BEAM PLASMA DYNAMICS 


The plasma beam that emerges from the thruster is a dense, charge neutral bea 
moving perpendicular to I. The physics of a plasma beam injected into a transverse 
magnetic field has been studied by many authors, both experimentally and theoreti 
cally (refs. 10,. LI, 12, 13). Also there have been authors who have considered th 
problem of ion thrttifrt in space, some specifically in the SPS context (refs. A, 

14, 15). 

Curtis and Grebowsky (ref. 14) argue that the beam polarises and £ - * J 

cancels the Lorentz torce. According to Curtis and Grebowsky (re . ) e ® 

density is always able to support the polarization field. In this case the beam 
simply passes out of the magnetosphere, depositing a thin non-propagatin^ sheath. 
While this is correct for a vacuum injection, when the field lines -nrea g 
wafer, shorted the plasma Is stopped (refs. 10 13). This Isexacty the wse U 

the magnetosphere. The beam dynamics paradigm of Chiu et al. (ref. . 

the barium release experiments (refs. 6, 12). In this picture the polarization 
field accelerates and polarizes the adjacent plasma, which in turn Pol" 1 *** ^ 
plasma adjacent to it. This electric field, which moves along the field line at 
the Alfven speed, transfers beam momentum to the ambient plasma and magnetosphere. 
SSen the AlfJIn wave reaches the ionosphere it drives dissipative Pedersen currents, 
and can be partially reflected (ref. 6). 

According to this model the plasma velocity decreases like e -t/x , where T l f 
the amount of time it takes for the Alfven wave to travel over as much ™* 8S un 

area as is causing the disturbance. Therefore, Chiu et al. (ref. 15) give that 


I dz P b /2V A P 0 
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where P Is the mass density (b refers to the beam, 0 to the ambient), V A is the 
Alfven speed, and the integral is along the field line. Calculations using realis- 
tic plasraaspheric and magnetic field models, give t « 10 seconds (ref. 16). Thus 
the beam can travel for distances £ 1000 km. 

Treumann et al. (ref. 17) have pointed out that as field aligned events 
out the polarization field, electrons cannot E * * drift across field lines to ne 


tralize the beam. They postulate that more currents parallel to tJ neutralize the 
beam head, generating a kinetic Alfven wave.- For the PORCUPINE Ion beam (of which 
more will he said), Treumann et al. (ref. 17) estimate electron drift velncities 
~V A » Due to these high drifts they argue that anomalous heating of electrons ener- 
gizes them to -20 eV to explain those hot electron s-Seen. in the PORCUPINE experiment 
(ref. 18). 

if we assume that the beam width is on the order of 10 km and the beam length 
~1000 km then the average neutralizing field-aligned currents (for the EOTV dis- 
cussed above) must be -16 uA/m 2 . For ambient electron densities of £ 10* m ’ this 
results in drift Speeds £ 100 km s”* . This is still half the thermal speed of a 0.1 
eV electron, so ion acoustic waves will be stable, but electrostatic ion-cyclotron 

waves, with w = a*, are unstable for V D > 5 * (ion thermal speed) (ref. 19) thus 
limiting the current. Therefore it is unlikely that the polarization field wil-l be 

completely shorted out and that the beam ions will b| charge-neutralized by a. combi- 
nation of both E x B drift of beam electrons across B and field-aligned currents. 

Haerendal and Sagdeev (ref. 8), writing on behalf of the PORCUPINE experi- 
menters, report on the injection of a 4A, 200 eV, Xe + , charge-neutralized plasma 
beam. This beam was injected at —72° to B in nine events ranging in altitude from 
196 km to 451 km. They report three stages of its beam expansion. The first is 
free expansion of the beam until the magnetic pressure starts to balance the dynamic 
pressure. The second phase is one of diffusive expansion, with the polarization 
field allowing for some motion across *> although the polarization field is rapidly 
shorted out by field-aligned currents and the beam is stopped. The third-phase is 
that of single particle motion. 


Considerable wave activity was also reported during the injection events (refs, 
18, 20, 21). Broadband ion-cyclotron harmonic waves were detected (refs. 18, 21) 
and it has been argued that the Drummond— Rosenbluth instability (ref. 19) is respon- 
sible (ref. 18). Given the much greater scale of the SPS injection it is reasonable 
to expect intense wave generation which could energize ambient particles. In par- 
ticular there is experimental evidence for the acceleration of thermal electrons by 
ion— cyclotron waves. Norris et al. (ref. 22) suggest that these waves in the mag- 
netosphere having f ) f|jg+ accelerate electrons with a clear bias parallel to q> 
the majority of the electrons heated to 20 eV. 

The field-aligned currents produced as a result of this are consistent with the 
view that such currents will play a central role in the dynamics of the beam, since 
these currents transfer momentum from the beam and drive Alfven waves. This is sup- 
ported by observations of considerably enhanced electron fluxes during the injection 
of a plasma transverse to B as reported by Alexandrov et al. (ref. 23). These 
authors also suggest that the observed magnetic disturbances are due to field- 
aligned currents and associated Alfven waves. 

The beam model presented below will deal with only the large-scale dynamics of 
the beam. It is assumed that the paradigm of Scholer (ref. 6) is essentially cor- 
rect, so beam velocity decreases as e~^ T . When the beam emerges from the thruster 
both nKT and (l/2)pv 2 are much greater than B ? -/2u 0 . Thus the beam will expand radi- 
ally outward from the beam axis, as if into a vacuum, forming a cone. This phase of 
the expansion continues until 


after which motion perpendicular to and ® is stopped, while expansion along the 
field- line continues essentially uninhibited (fig. 1). The beam travels In this 
manner until 

i pv 2 a tt?/2u 0 (2-3) 

at which point the geomagnetic field becomes the dominant influence in the argon 
plasma's motion. The argon's subsequent motion can -then be -followed 4>y the adia- 
batic theory. 


In the initial phase of the plasma beam we are dealing with a vacuum expansion. 
It is assumed the beam Is charge-neutral and colllsionless. In the very early his- 
tory of the beam it is certainly collisional, the plasma rapidly thermallzlng, but 
by ~100 m downstream the mean^free path is of the order of the beam size. We also 
neglect momentum loss to Alfven waves, assuming t << t during the first (vacuum) 
phase of the beam expansion. 

To represent the plasma that emerges from the thruster we write 


f(X, Y, V 2 , t = 0) 


litKT 5 ^ V z ” V B^ exp 
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where % = Xi + Yj 

+ V = V x i + Vyj 
Vjj = beam velocity = Vj$z 
R = 1/3 r 0 

Liemohn et al. (ref. 4) give KT - 10 eV, which for a 2 keV beam gives an effective 
beam divergence angle ~$°. This same 8° spreading angle is reported by Cybulski et 
al. (ref. 24) in flight tests of the SERT I ion thruster. To confine the plasma to 
the thruster at t » 0, R is set to one-third the thruster-radius, rjj. 


We can then write the colllsionless, vacuum Boltzman equation 
|| + $ • Vf =* 0 

The solution to equation (2-5) is readily found to be 


(2-5) 


f(X, V, v z , t) 



(X ~ frt) 2 i 
2R 2 J 
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The constant, n 0 , is given by the normalization condition 

/ fd 2 v dx V t = N = !£• 

B e 


(2-7) 


where the integral along the beam axis is replaced by multiplication by Vgt. From 
equation (2-7) one finds 
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The density of ions is given by 
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2R 2 (ra/KT) + ( t 2 / R 2 ) 


which for t X )0~ 4 8 y lelds » u9tn 8 equation-(2-8) 
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The end of the vacuum expansion phase is given by equation (2-2) using 
'2-10) evaluated at the edge of the beam, which is r « V Th t. The time at which th 
condition is satisfied for the EOTV iti question in a dipole field is 


8.83 x 10” 3 L 3 sec 


(2-11) 


where L is the magnetic shell parameter. After time t 0 the beam continues Jo spread 
along $ and so the density decreases like 1/t. But the velocity is also going down 
^t/T and 80 the density must go like e^ x to conserve particles. So we may 


as e 


write the central beam density for t > tg as 

t n(r = 0, t = t Q ) _ t/l 
n(t > t n ) 1 6 


( 2 - 12 ) 


Using this density in equation- (2-3) yields an equation for the time of transi- 
tion. from beam motion to adiabatic motion, which is. 

L 3 = 0.216 t e t ^ T (2-13) 

where L is the dipole shell parameter. The fraction of energy the beam ions retain 
is then e- 2t '\ the rest of the energy being transferred to the magnetosphere and 
ionosphere. As stated before, t - 10 sec in the plasmasphere (ref. 15). Figure 
gives the energy loss as a function of L for x - 10 sec from L * 2.5 to L - 3.5. 

To estimate the energy density deposited in the ionosphere by the beam we find 
that at L - 2 roughly 3/4 of the beam power is lost to Aifven waves. From equatio 
(2-11) we find that the beam width ~ few km, and the beam length < 1000 km. This 
eives an area of < 10 4 km 2 , which mapped down to the ionosphere (with a dipo 
field^resulta in~ionospheric power densities of -lOT* W/m 2 (assuming the bulk of 
the energy is deposited in the ionosphere). This is considerably larger r than ; auro- 
ral nower densities of~10“ 2 W/m 2 (ref- 25), and even if we assume only 10% of the 
power is absorbed by the ionosphere there would still be power densities equivalent 

to aurorae. 


CONVECTION OF THE ARGON PLASMA 


Once the argon plasma’s transition from beam to individual particle motion is 
accomplished the subsequent motion is determined by the local magnetic and 
fields. This motion is most easily followed using the guiding center approximation. 
The tlo first orde?drift velocities are the I * 1 and the gradient/curvature drift 
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velocities* The motion of the plasma la gi ven by 


dx 

dt 


m 


1 X ft c 1/2) 'mv 2 
B?- B 


(l + cos 2 «) 


IjLUlL 

eB 2 - 


(3-D 


Ignoring time dependence, ft and $ are functions of & and this equation must be 
solved numerically* 

The coordinate system Is as follows: x la the antlsunward direction, y along 

the dawn meridian* For a magnetic field the model of Mead (ref* 26) is used with 
the magnetopause set at r « 10 R$. In the equatorial plane, B just has a z compo- 
nent: 


B *» B 

z 


3*11 x 1(T 5 
R 3 


+ 2.515 x 1(T 8 - 2.104 x 10" 9 R cos 6 


(3-2) 


where R the earth radii,. 6 is the local time and B in tesla. 


There are three components to the electic field: the convection field, the 

corotation field and the self-electric field of the argon plasma, which will be dis- 
cussed in more detail below. The convection field is approximated by a constant, 
dawn to dusk, 0.2 mV/m electric field. This corresponds to a ■ 50 kV cross-polar 
potential drop mapped out onto, a 40 R^ magnetosphere. The corotation field is given 
by 


-V -r- 

E = -V x B 
cor 



x r 


x B Dipole 


(3-3) 


The above model has some obvious shortcomings. The- greatest of these is the 
assumption of a uniform convection field. In addition, there is an inconsistency in 
making the approximation in equation (3-3), then using it in the drift equation, 
since the expression for $ has non-dipole components. The proper way to calculate 
the corotating field is to calculate v x B for the field line in question. This was 
not done due to the limitations of the Tektronlcs 4052 minicomputer which was used. 
However, this rough model should give a somewhat reasonable approximation to the 
general features of the convection of injected argon. 


Chiu et al. (ref. 15) give the time fraction spent, from LEO to GEO in figure 3. 
Using this with the above model one finds that —3/4 of the injected argon is trapped 
in the plasmasphere. However, a more complete picture of the argon convection 
requires the inclusion of the electric field generated by the argon filled flux tube 
itself. The gradient drift current in the flux tube (or plasma blob) gives rise to 
Birkeland currents, which close in the ionosphere, if neighboring flux tubes cannot 
satisfy continuity of current. This current system is illustrated in figure 4. 


Assuming that equal amounts of current go to the northern and southern hemi- 
spheres current continuity gives 


— V • 5 
2 J leq 
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where is the current/ length in the equatorial plane, J|| e q is the current/unit 
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area out of the equatorial- plane* and J|ii on . the current/unit area- into the iono- 
sphere. This current closes la the ionosphere, in which case we can take the 
divergence of Ohm's Law and, using equation (3-4), ..Obtain 

B 

?*<S-VV>-ivj (jiffi) (3-5) 

H eq 

♦ •> 

Assuming J ton to be 1 B, and considering only the Pedersen conductivity, the 
above becomes 


1 ^ I). 

= 2 " ^ * ^ ea (g °°) » V. ® horizontal 7 in ionosphere (3-6) 

eq 

This approach, developed in part by Vasyliunas (ref. 27) and Wolf (ref. 28), allows 
the calculation of the ionospheric potential set up by the flux tube. This poten- 
tial can be then mapped out along field lines (assuming they are equlpotentials) to 
give the potential, and electric field. In the equatorial plane. 

We assume the argon density to be constant throughout a flux tube of radius a, 
and zero outside of the flux tube. It is also assumed that there is a uniform, 
background current density. For this case, in radial coordinates centered on the 
flux tube-, we find 

V • J ieq = K cos 9 6(r - a) (3-7) 

where K is a constant. As a further simplification we assume we are dealing with a 
circle in a locally flat ionosphere, so equation (3-6) becomes 

V * V = f“ cos <|» 6(p - a') (3-8) 

P 


where (p,^>) are the ionospheric coordinates whose origin is the field line threading 
the center of the plasma blob, and K' , a' are constants. The solution is obtained 
ia a straightforward fashion to yield, in the ionosphere 

v ■ fro - cos * » P > a ' 

P 


K'p . 

-Jjp cos i|i, 

P 


P < a' 


(3-9) 


Mapping the resultant electric field out 
field lines yields 

t K.1L. 

Blob 2ZpX 3 / 2 


to the equatorial plane along dipole 


(3-10) 


where Is local time and X is the L shell of the^cen^er of the argon flux tube* 

This electric field, for negative K', results in Ex B motion radially outwards from 
earth. We also note that the field outside the blob is that of a dipole. 


The constant K' is related, by equations (3-6) and (3-7) to the divergence of 


the net current, Which Is 

7 • V - ( I J bcL- l^cc 11 cos 6 S(t - a) (3 - ll) 

where |Jq(;| Is the-graiHeitf-cervature current, and |J,gj ta- Che background current 
present in neighboring flux tubes* 

We can consider two extreme cases concerning the condition of the argon in 
the flux tube: very strong pitch angle scattering so that the distribution is iso- 

tropic, or very weak scattering so the plasma mirrors at the equator (~ the state of 
injection). In the later case 
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where N is the nuraber/area and p is the magnetic moment, to find P one can use the 
beam model of Chapter 2 to obtain the ion's perpendicular energy when they Start 
convective motion. For a dipole field this gives 




B 


(3-13) 


In the former case the flux tube can be treated as an ideal gas »_ in which case 
an adiabatic energy invariant X can be defined- (ref . 29) 


e k - - » ( / */»)-*/• 

One can— also define a number invariant (ref. 29) 

n = n / ds/B * N/B = ///flux; n = ///m 3 
then the gradient-curvature drift current is given by (ref. 29) 
J GC = nXz X 7 eq (/ ds/B)-2/ 3 

Using equations (3-14) and (3-15) the above can be written 
ez x 7 (/ ds/B) 


GC 


(3-10 


o-m 


C 3-lo > 
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where e is the energy density of the flux tube. This form is convenient since 
Williams (ref. 30) gives the quiet time energy density to be ~10” 3 Joules. 'm- 
from L ~ 2-5.5, thus we can calculate the background gradient-curvature drift cur- 
rent, 5g G . 

To proceed we need to calculate n, X, /ds/B, and V(/ds/B). The flux tube vol- 
ume is given by G.-H. Voigt (personal communication, 1982) to be 


/ ds/B = - LI 1 -— [l + ^ L + ® L 2 + L 3 ] 


y 


(3-lS> 


the gradient of /ds/B can be obtained by straightforward differentiation of equation 


683 


(3-18). In. addition, X may be obtained by using equation (3-i4) with the ion H no 
tto energy given by the beam model In Chapter 2. ' ' 1 C C 1 kl 

number^denalty^ 1^^’ 'if fr ° m aquat4on (3 ~ 15 ) if ^ the area 

i ucuniry, is Known, u i i» the Ion- beam current of. the thruatom m/» v 

la the orbital velocity of the vehicle then N la given by E0T ^ 

N • i ' , ( qV kotv <I t ) (3-19) 
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the bulk of the ti£n *TtZ\T oTlZ JT™ ? the J COOVeCtlon ° odel is that 
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10-100 (ref 93 n . mailer than the ambient density by a factor of 
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phenomena In the tall. an ws c hanglng the time sealea for dynamic 
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radialLy outward they produce an inertial drift current oppositely directed to the 
gradient-curvature drift current. So the flux tubes will accelerate at the rate 
needed to cancel the gradient-curvature current, thus Satisfying continuity, until 
the velocity of the Ar*" has- reached is that of the E * n velocity of the self-, lec- 
tric field, at which point it will- not need— to- a cc o lor -o t e to- close- the current. 


The inertial dd£c current/ length la 


5 * - Nroft x B 
I eB 2 


0 - 21 ) 


where N is the equatorial area number density and m Is the argon mass. Equating 
this to the gradient drift current and solving for 3, the acceleration, gives (for 
B " B„/L 3 ) 


& S 


3UB n 
0 * 

mi^L 4 


(3-22) 


For argon injected + at L « 2 equation- (3-22) yields a ~ 5 km/s 2 so within a short 
time St ~ Eg^oB * &/b 2 . Thus, except for the first moments of the argon drift. 
Inertial currents need not be self-conslstently included. Gravitational and centri- 
fugal drifts may also be neglected. 


CONCLUSION 


We have seen that the operation of powerful ion thrusters In the SPS context 
(or for that matter in any Space Industrial project of such scale) is expected to 
have a wide variety of effects. Two distinct plasmas are injected: a 2 keV beam 

plasma and a thermal charge-exchange plasma. The thermal plasma will be a signifi- 
cant addition to the heavy ion content of the plasmasphere. Immediate effects of 
the beam injection include wideband ion-cyclotron wave generation and field-aligned 
heating cc electrons. Alfven waves will transfer beam momentum to the ambient 
plasma and ionosphere. In the ionosphere, power densities will be of the order of, 
or bigger than natural auroral power densities and so the EOTV will in effect create 
an artificial aurora on the order of 10 2 km long. These induced aurorae will pro- 
bably have localized, disruptive effects on communications, and could also affect 
pdwer transmission lines (ref. 32). 

Once the beam Ions start to convect, they move radially outward until the argon 
flux tube energy density Is equal to the background energy density. Therefore most 
of the argon will convect out of- the plasmasphere. The field-aligned currents gen- 
erated could cause further ionospheric disturbances. The argon will then enter the 
general convective cycle and be distributed throughout much of the magnetosphere. 
This significant number of heavy Ions could substantially alter dynamic quantities 
like the Alfven speed. The dispersion relation for waves would also be altered. 

New cutoffs and resonances should appear along with new wave-particle phenomena. In 
fact, Chiu et al. (ref. 13) have shown that the presence of Ar + damps the Instabi- 
lity which precipitates MeV electrons, therefore those levels should rise. 

Thus, we see that the operation of an argon-ion propelled orbital transfer sys- 
tem will substantially modify the magnetosphere-ionosphere system. These raodlflca- 



tionn, while being substantial, do not seem to be, In and of themselves, so drastic 
as to rule out this transportation system, historically, whenever man enters a new 
environment he modifies hla culture, technology, and himself, while in turn altering 
the environment, both by his activities and to suit his needs. The magnetosphere in 
no exception, and as man becomes more 1 ,; wo J wd - in - sp ace In the hear future, he will 
begin to modify it.- 
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(Adapted from Freeman et al 
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